Noise-induced Synchronization in Small World Network of Phase Oscillators 



Reihaneh Kouhi EsfahaniQ Farhad ShahbaziQ and Keivan Aghababaei Samanl§ 
Department of Physics, Isfahan University of Technology, Isfahan 84156-83111, Iran 



A small-world network (SW) of similar phase oscillators, interacting according to the Kuramoto 
model is studied numerically. It is shown that deterministic Kuramoto dynamics on the SW networks 
has various stable stationary states. This can be attributed to the defect patterns in a SW network 
which is inherited to it from deformation of helical patterns in its parent regular one. Turning on 
an uncorrelated random force, causes the vanishing of the defect patterns, hence increasing the 
synchronization among oscillators for intermediate noise intensities. This phenomenon which is 
called stochastic synchronization generally observed in some natural networks like brain neuronal 
network. 

PACS numbers: 05.45.Xt 87.19.1c 89.75.Hc 



INTRODUCTION 



Noise is usually considered as a source of disturbance 
against the main signals in laboratory as well as natural 
systems. Nevertheless, the interplay between the ran- 
domness, created by the noise, and non-linearities may 
lead to the enhancement of regular behaviour in some dy- 
namical systems Stochastic resonance 0], coherence 
resonance 13], noise-induced transport noise- induced 
transition [5j and noise-induced collective firing in ex- 
citable media Q are examples of such a novel phenom- 
ena. Being noisy, nature takes the advantages of these 
machanisms in order to employ the random fluctuation 
as agent of self-organization. This is the main reason 
that why living systems work so reliable in spite of the 
presence of various sources of noise. Brain neurons are 
examples of biological systems in which the source of 
random fluctuations is the background synaptic noises 
caused by highly fluctuating inputs coming from thou- 
sands of other neurons connected to a given neuron Q. 
However, this noise plays a constructive role in regular 
spiking of the individual neurons and also increasing the 
synchronization among the clusters of connecting neu- 
rons [8] . Synchronous spiking among a subset of neurons 
plays an important role in more efficient propagation of 
activities from a group of neurons to another [9(. Fur- 
thermore, there are also some controversial idea on en- 
coding of information about stimuli thorough synchrony 
in oscillatory activity of neurons flpj ] . Another phenom- 
ena in which noise-induced synchronization takes place, 
is gene regulatory processes in systems such as quorum- 
sensing bacteria, in which noise originates from the small 
number of molecules involved in the related biochemi- 
cal reactions [ll|. Collective dynamical behaviours, like 
synchronization, can be simulated by systems of coupled 
non-linear oscillators. One of such models has been pro- 
posed by Kuramoto, which consists of a set of oscillators 
with fixed amplitude (phase oscillators) mutually coupled 
with sine of their phases difference [12|. The stochastic 



Kuramoto model has been studied on the globally con- 
nected [lU and also on scale-free (SF) and Erdos-Renyi 
(ER) random networks [13] • Analytical results on the all- 
to-all network show that for a given distribution of intrin- 
sic frequencies of oscillators, a minimum value of coupling 
is needed for them to become synchronized. Deriving 
the synchronized system by an uncorrelated white noise, 
causes the synchrony between oscillators declines mono- 
tonically by increasing the noise strength. The same re- 
sults have been found on numerical integrations of the 
stochastic Kuramoto model on the ER and SF networks. 
The difference is that the synchronized state in SF net- 
works persists more against applying the noise with re- 
spect to the ER and all-to-all networks flij. In an effort 
to search for the reason of synchronization between crick- 
ets, Watts and Strogatz found out that many systems in 
nature possess the properties of the small-world (SW) 
networks 15|, [16| . Short mean path between the nodes 
and high degree of clustering are the two main features 
of SW networks. Former is a characteristic of random, 
while the latter is feature of regular networks. It has 
been found that the presence of random short-cuts, may 
lead to noise driven ordering phenomena such as stochas- 
tic resonance [13] and coherence resonance 18] in SW 
networks. Motivated by recent discoveries revealing SW 
topology of brain neural networks [l9| and also noise in- 
duced regulatory behaviours in such a networks [8|, we 
study the effect of random force on dynamics of SW 
network of a set of similar phase oscillators coupled to 
each other based on Kuramoto model. We will show 
that in this system, for intermediate noise strength, the 
synchronization among the oscillators is increased. The 
rest of the paper is organized as follows. In Sec. II, 
we present the results of numerical integration of deter- 
ministic Kuromato model on regular and SW networks. 
Investigation of Stochastic Kuramoto model driven by 
uncorrelated white noise is done in Sec. Ill and Sec. IV 
is devoted to summary and concluding remarks. 
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KURAMOTO MODEL ON COMPLEX 
NETWORKS 

In this section we introduce the Kuramoto model and 
numerically investigate its steady state solutions on ER, 
SF and SW networks. Consider a set of phase oscillators, 
residing on the top of the nodes of a network. Their 
phases and intrinsic oscillation frequencies are given by 9i 
and Ui, respectively. According to the Kuramoto model 
the dynamics of these phase oscillators is given by the 
following set of coupled differential equations: 
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where K is the coupling strength, N is the number of 
nodes and a,-j is the element of adjacency matrix (etjj = 1 
if nodes i and j are connected and a,j = otherwise). 

The synchronization of Kuramoto model on SW net- 
works, for random distribution of Wj has already been 
studied by Hong, et al [20| . They showed that small frac- 
tion of shortcuts is enough for both phase and frequency 
synchronizations on SW networks, in spite of absence 
of any synchronization on regular ones. In our work, 
we assume that all the intrinsic frequencies are the same 
(uji — wo) i therefore moving to a reference frame in which 
ujq = 0, simplifies Eq.([T]) to: 
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For comparison of the solutions of Kuramoto model in 
these three types of network, we need to construct them 
with equal number of nodes and edges. For building a 
SF network with average connectivity (k) = 2m, we use 
the Barabasi- Albert algorithm [21] . Starting from mo 
initial connected nodes, one attaches a newly entering 
node to m < mo elder ones with probability propor- 
tional to the degree of the present nodes. An ER ran- 
dom network with N nodes and the same average degree 
per node ((k) = 2m), is simply produced by connect- 
ing randomly chosen pair of nodes with Nm edges [22j |. 
To construct the SW network, we use Watts-Strogatz 
(WS) algorithm 15|. Starting from a regular network 
with N nodes and k — 2m edges for each node, we 
rewire each edge randomly with probability p. Choos- 
ing 0.005 < p < 0.05, this process converts the initial 
regular lattice to a complex network with a small mean 
path length and large clustering coefficient, character- 
istics of SW networks. Starting from a randomly dis- 
tributed initial phases 0,(0) (which is selected from a box 
distribution in the interval [— tt,7t[), The set of coupled 
differential Eqs.@ are integrated from t = to a given 
time t with the time step dt, using Euler method. This 
method enable us to compute 9i(t) and to determine the 
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FIG. 1: (Colour on-line) Order parameter (r) versus time for 
SW, SF and ER networks with N = 1000 and (k) = 10. 
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FIG. 2: (Colour on-line) Time dependence of the order pa- 
rameter (r) for a SW network with iV = 1000 and (k) = 10. 
Different curves denote different initial conditions. 



synchrony among the oscillators at any time, we define 
the following complex order parameter: 



1 N 

-Y 



(3) 



where < r(t) < 1 indicates the degree of synchroniza- 
tion in the network and ip is the phase of the order pa- 
rameter. 

Fig.([T]) shows the temporal variations of the r(t) on 
the three types of network with N = 1000 and (k) = 10. 
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of Eqs.Q have to satisfy the following conditions 231 ] : 
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FIG. 3: (Colour on-line) steady state phase configurations 
in 4 helical patterns of regular network and corresponding 
steady states of SW networks with N — 1000 and (k) — 10. 
(b) A = 1000, (c)A = 250, (d)A = 100 and (e)A = 50. A is the 
wavelength of helical states in regular network. 



To obtain these plots, time step is set to dt = 0.01 and 
averaging is done over 100 realizations of initial phase dis- 
tributions for a fixed network of each type. The rewiring 
probability for constructing the SW network out of reg- 
ular one is chosen to be p — 0.04. As can be seen, the 
oscillators on ER and SF networks immediately reach to 
a fully synchronized state (r = 1), while in the case of 
SW network they more slowly go toward a partially syn- 
chronized state with r « 0.7. These results show that in 
contrast to ER and SF networks, the structure of steady 
states of the Kuramoto model on SW networks are more 
complex. To get more insight on the nature of steady so- 
lutions, we check the dependence of order parameter on 
the initial phase distributions. Fig. shows temporal 
variation of r on a SW network for five different initial 
conditions. As can be seen from this figure, the steady 
state is highly dependent on the initial phase distribution 
in such a way the r(oo) reaches several values between 
and 1. In what follows, we discuss that the sensibility of 
dynamics to initial conditions is indeed inherited to SW 
networks from their regular network parents. 

It is easy to show that the stable stationary solutions 
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provided that the phase difference between any two adja- 
cent oscillators be less than tt/2 (i.e,A9ij = 9i —9j < ir/2 
if ciij = 1). These solutions can be put in two categories 

i) Full synchronized state with r = 1 (A6*y = for any 

ii) Phase-locked state with regular arrangement of phases 
around the phase circle with non-zero phase difference 
A9, for which r = 0. 

The phase-locked states represent helical wave phase 
modulations and their number depends on N and k. 
For instants, in the case of A = 1000 and k = 
10, there are 10 of such states with nearest neigh- 
bour phase differences A9" n = 27r/A Q , in which A" = 
20,25,40,50,100,125,200,250,500,1000, are the wave- 
lengths of the helical states indicated by a = 1,2, •••10, 
respectively. 

The stationary phase configuration of all nodes, cor- 
responding to the initial conditions in Fig. ([2]), are plot- 
ted in Fig. ©, both for the regular and its offspring 
SW networks. This plots corresponds to the helical pat- 
terns with phase differences A = 1000, 250, 100, 50, de- 
noted in Fig.© by indices, b, c, d and e, respectively. 
This figure shows, rewiring a regular network with phase- 
locked state, deforms its helical pattern to an inhomoge- 
neous state in the subsequent SW one. Therefore, a SW 
network possesses various stable stationary states whose 
number equals the number of helical patterns in its par- 
ent regular network. 

The local structure of the steady state can be better 
clarified by the correlation matrix D defined as 24[: 
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cos(6i{t) -6j(t))dt, (5) 



in which t r is the time needed for reaching to stationary 
state. The matrix element — 1 < D^j < 1 is a measure of 
coherency between each pair of nodes. In the case of full 
synchrony between i and j (9i = 9j) the correlation ma- 
trix element is Dij = 1 and in the case of anti-phase lock- 
ing (9i — 9j = 7r) , the value of matrix element is = — 1 . 
Fig. (0| represents the density plot of correlation matrix 
elements for the four steady states of regular and SW 
networks corresponding to Figs. ([2]) and ((SJ. This plots 
clearly show the inhomogeneous structure of the helical 
patterns before and after rewiring of the regular network. 
The correlation matrix represents strip structures in its 
density plot for helical states in regular network and the 
width of the strips are proportional to the wavelength of 
the helices. One can also observe from this plots that 
converting the regular network to SW, the helical pat- 
terns are substantially affected, provided A being large. 
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FIG. 4: (Colour on-line) Density plot of correlation matrix 
elements (Dij) for 4 helical states of regular network and cor- 
responding stationary states in SW network with N = 1000 
and (k) = 10. (b)A = 1000, (c)A = 250, (d)A = 100 and 
(e) A = 50. A is the wavelength of helical states in regular 
network. 



The strip structure of matrix D is almost preserved for 
small wavelengths, indicating the small wavelength he- 
lical patterns, despite of little deformations, are stable 
against rewiring of network. For large wavelength pat- 
terns of regular network, the majority of nodes in the 
corresponding SW phase configurations are synchronized 
with each other, however there are some isolated nodes in 
Anti-phase locking with the rest. These isolated nodes 
are topological defects and induce spiral phase textures 
around them, in such a way that the phase of surround- 
ing oscillators varies continuously from to 7r, by getting 
away from these nodes. The number of these defects in- 
creases by decreasing the wavelength of corresponding 
helical pattern. For example, it can be seen from Fig.© 
that for A = 1000 there is one while for A = 250 there 
are four point defects. Once the structure of the steady 
states of deterministic Kuramoto model on SW network 
is known, it would be interesting to investigate the effect 
of noise on such states. 
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FIG. 5: (Colour on-line) Order parameter versus noise in- 
tensity for SW, SF and random networks. Number of nodes 
and mean degree for the three networks are N = 10000 and 
(fc) = 10, respectively. 



EFFECT OF RANDOM FORCE 

A network of oscillators could be plagued by some ex- 
ternal random forces. The effect of such forces may be 
modelled by an uncorrelated white noise (r]i{t)) applying 
to all nodes. Adding this noise to Eq. ([2]), we have: 
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where (^(i)) = 0, (Vi{i)Vj(t')) = 2D8{t-t')8 ij with D 
being the variance or intensity of the noise. In our nu- 
merical work, we choose a box distribution in the inter- 
val [—w/2,w/2] for r], so that its variance is equal to 
D = w 2 /24. It can be shown that by proper rescaling of 
the time variable, the effect of parameters D and K can 
be included in a single parameter g 2 
ing the dynamical equations to: 
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where r = Kt is the rescaled time variable and £i(r) := 
r)i(t)/g is a random variable in the interval [—1/2, 1/2]. 

The numerical integration of Eq.© is carried out by 
employing Euler method for its deterministic part and 
Ito's algorithm 25] for the stochastic part. Fig.© rep- 
resents the variations of stationary order parameters 
(r(oo)) versus rescaled noise intensity g, for the three net- 
work types SF, ER and SW. In obtaining these graphs, 
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FIG. 6: (Colour on-line) Density plot of correlation matrix 
elements (Dij) for a steady state with four point defects on 
SW network and different noise intensities, (a) g = 0, (b) 
g = 2, (c) g = 3, (d) g = 4, (e) g = 5, (f) g = 6, (e) g = 7 and 
(f) g = 8. Number of nodes and mean degree for the three 
networks are N = 10000 and (k) = 10, respectively 



FIG. 7: (Colour on-line) Probability distribution function of 
correlation matrix elements (Dij) for different noise intensi- 
ties: (a) g = 0,(b)g = 2, (c) g = 3,(d)g = 4, (e) g = 5, (f) 
g — 6, (g) 3 = 7 and (i) g = 8. Number of nodes and mean 
degree for the three networks are N = 10000 and (k) = 10, 
respectively. 



we retained the networks fixed and averaged over 100 
different realizations of noise. 

By inspecting this figure one can extract two essen- 
tial results: (i) While The order parameter begins from 
a value less than 1 for SW network its variation is much 
smaller by increasing noise intensity, relative to ER and 
SF networks. Hence , at small noise intensities, the par- 
tially synchronized state in SW network is more robust 
against applying the noise than fully synchronized ones 
in the ER and SF networks. However, at large noise in- 
tensities the situation is vice versa. The critical coupling 
(g c ) at which the synchrony disappears among the os- 
cillators is the greatest for SF and the smallest for SW 
network. Therefore the coherent state in the SF lasts 
longer against noise than ER and SW with the same av- 
erage degree and coupling constant. The persistence of 
synchronized state in SF network has been argued to be 



related to the existence of few nodes with very large num- 
ber of connections (hubs) in this type of networks 14j . 

(ii) Increasing the noise intensity, the coherency among 
the population of phase oscillators destroys monotoni- 
cally in SF and ER networks. Nevertheless, in SW net- 
work, there is an interval of reduced noise intensity i.e 
5 < g < 7, within which the synchronization enhances 
by increasing noise strength. Therefore in this region, 
instead of destroying the synchronization, noise promote 
the synchrony among oscillators. 

The noise-induced synchronization is also called 
stochastic synchronization and its occurrence in SW net- 
works can be explained in terms of defect patterns in the 
steady states of Kuramoto model. Fig. © represents the 
evolution of correlation matrix density plots versus re- 
duced noise intensity, g, for a specific steady state of SW 
network with four topological defects. It can be seen in 
this figure that turning the noise on, the defects resist 
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against the noise up to g ~ 4 and for g > 4 they begin to 
disappear until g ~ 6 where they vanish completely. Dis- 
appearance of defects enhances the homogeneity in the 
system and so the synchrony among the oscillators. This 
is more apparent in probability distribution of correla- 
tion matrix elements (p(Djj)) shown in Fig. ([7]). As can 
be seen in this figure, p(D) has two peaks at D = 1,-1 
for g — 0. Increasing the noise intensity the two peaks 
move toward each other and at the onset of stochastic 
synchronization, g ~ 6, they emerge in one peak. At this 
point the variance of p(Dij ) reaches to its minimum and 
again rises by increasing the noise strength. 

CONCLUSION 

In summary, we found that a SW network of similar 
phase oscillators communicating with each other by Ku- 
ramoto coupling shows novel behaviours. Unlike ER and 
SF networks, this system fails to reach full synchronized 
state. Moreover, driving it by an uncorrelated with noise, 
reveals the occurrence of stochastic synchronization, a 
phenomenon through which a random force induces syn- 
chrony among the oscillators. We discussed that the rea- 
son for this phenomenon is lied in the stable helical pat- 
terns in the regular networks from which the SW ones 
is built. Rewiring of a regular network of similar phase 
oscillators with periodic helical pattern, ends to complex 
inhomogeneous states in the resulting SW network. The 
existence of such stable inhomogeneous patterns in SW 
network, appearing some times as topological point de- 
fects and also as aperiodic helical patterns, prevents the 
network from reaching to full synchrony. These patterns 
persist against applying the noise for small noise intensi- 
ties. However the external random forces with moderate 
strengths are able to destroy these patterns in favour of 
more homogeneous states, hence enhance the synchro- 
nization among oscillators. At the end, we hope that 
our finding that such a simple model representing such 
a novel behaviours may shed light on the fact that why 
SW networks are so ubiquitous in natural systems. 
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